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By 

C, C. Chamis* and M- D. Minich ** 

Lewis Research Center 
National Aeronautics and 
Space Administration 
Cleveland, Ohio 


Abstract 


A theoretical investigation was performed to determine the structural re- 
sponse of a fiber composite airfoil typical of those encountered in high- 
tip speed compressor fan blades when subjected to load conditions antici- 
pated in such applications. The analysis method used consists of composite 
mechanics embedded in pre- and post -processors and coupled with NAS TRAN. 

The load conditions examined include thermal due to aerodynamic heating, 
pressure due to aerodynamic forces, centrifugal, and combinations of these 
loads. The various responses investigated include root reactions due to 
various load conditions, average composite and ply stresses, ply delamina- 
tions, and the fundamental modes and the corresponding reactions. The 
results show that the thermal and pressure stresses are negligible compared 
to those caused by the centrifugal forces- Also, the core-shell concept 
for composite blades is an inefficient design (core plies not highly stressed) 
and appears to be sensitive to interply delaminations. The results are 
presented in graphical and tabtilar forms to illustrate the types and amount 
of data required for an analysis such as this, and to provide quantitative 
data associated with the various responses which can be helpful in design- 
ing composite blades. 


*Aerospace Engineer, Materials and Structures Division 

**Associate Professor, Department of Civil Engineering, The Cleveland State 
University, Cleveland, Ohio 
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INTRODUCTION 

The recent national awareness of the need to conserve our natural re- 
sources provides a strong motivation for efficient material utilization 
and the development of new design concepts for structural components in 
general. One particular area where efficient material utilization and new 
design concepts could help conserve natural resources is in the use of fiber 
composites in high-tip-speed airfoils for compressor fan blades in aircraft 
engines. Increasing the airfoil tip speed permits reduction of the number 
of stages in the compressor to attain a given thrust. By using advanced 
composites, higher airfoil tip speeds and lighter weights can be attained 
compared to the use of state-of-the-art materials. Both higher tip speed 
and lighter weight airfoils have an accumulative effect on reducing the 
amount of material used in an aircraft engine. 

Advanced fiber composites application to high-tip-speed airfoils is 
currently under investigation by the NASA Lewis Research Center under con- 
tract NAS3-15335 with Pratt and Whitney Aircraft, East Hartford, Conneticut. 
Another application is reported in reference 1. 

A computerized analysis capability which can be used to analyze fiber 
composite airfoils is described in reference 2, along with analysis results 
for a specific design. However, the structural response of composite airfoils 
due to various loading conditions has not been investigated. 

Determining the structural response of the airfoil due to various load- 
ing conditions is useful in several respects. For example: (1) The airfoil 
stiffness and thermal properties may be used to obtain an a priori assessment 
of the deformation and stresses in the airfoil when subjected to various load 
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conditions- (2) The stresses, deflections, and root reactions produced by 
thermal, aerodynamic pressure and centrifugal loading conditions can be exam- 
ined to determine which is the dominant loading condition. If a dominant 
loading condition is established, then only this loading condition need be 
used in preliminary designs, thereby saving considerable man-hours of effort 
and computer time. (3) Tip deflections due to combined loading may be studied 
in order to determine the need for offsetting unacceptable radial growth, 
uncambering and untwisting- (4) End reactions due to vibrations may be studied 
to determine whether high cycle fatigue may be a design problem at the root 
at t achment . 

The primary objective of this investigation is to use the computerized 
capability (ref. 2) to theoretically examine the structural response of a com- 
posite airfoil, typical of those encountered in high-tip-speed compressor fan 
blade applications. The airfoil is assumed to be subjected to the following 
loading conditions: thermal, dynamic pressure, centrifugal, a combination of 
these, and vibration. The structural response of interest includes: airfoil 
stiffness and thermal properties, displacements, composite and ply stresses, 
root reactions, thro ugh- the -thickness shear forces and stresses, interply de- 
lamination, airfoil fundamental frequencies and their respective root reactions. 
The centrifugal load stiffening effects and vibration mode shapes are exten- 
sively discussed in reference 2 and, therefore, are considered in this inves- 
tigation only briefly. A secondary objective is to provide quantitative 
guide data for designing airfoils for high-tip-speed compressor blade appli- 
cations . 

The airfoil geometry, aerodynamic loads, and rotational speed selected 
are representative of the preliminary design of a high-tip-speed blade, and 
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were obtained from the contract previously mentioned • The laminate config- 
uration is of the core-shell type: 0® plies for the core and +45® plies for 
the shell. Briefly this capability consists of composite micromechanics , 
macromechanics * and laminate analysis embedded in pre- and post-processors 
and coupled with NASTRAN. 

INPUT DATA 


The input data for the computerized analysis consists of the following 
items : 

1. Airfoil geometry defined at several cross sections 
along the airfoil span 

2. Airfoil operating temperature 

3. Airfoil dynamic pressure difference between pressure and 
suction surfaces 

4. Composite system 

5. Airfoil laminate configuration, i.e. number of plies 
and ply orientation 

6. Finite element representation including load case selection, 
type of solution (static, dynamic) and boundary conditions* 

For the specific blade under consideration, the airfoil geometry was 
specified as input at three cross sections; root, midspan and tip. The air- 
foil operating temperature and pressure difference were specified as input at 
the same sections for convenience. Pertinent data, obtained from the contract 
previously mentioned, are summarized in Table 1 which is a computer output of 


the input data. 
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The composite system selected was Modmor II (Mod-II) graphite/PMR poly- 
mide. The ply thickness used was 0.007 inches at 0.57 fiber volume ratio 
and zero voids. The laminate configuration for Che airfoil was of the core^ 
shell version, symmetric about the midplane. The ply orientation was as 
follows : 60 percent of the plies at the airfoil maximum thickness were ori- 
ented at 0® (core plies) and the remaining 40 percent were at +45° (shell 
plies). The 0° direction is parallel to the span or radial direction. These 
input data are summarized in Table 2 which is a computer output of the com- 
posite type and ply geometry. In Table 2 the ply starting (%I) and stopping 
(%F) points (ply contours) are given as distance percentages along the co- 
ordinate axes of the airfoil. See also Nomenclature. Note the type of for- 
mat in Table 2 was selected to accomodate usage of hybrids in both the same 
ply level and through the thickness. This type of format also accomodates 
specification of positioning of leading edge devices. 

The last part of the data input is the finite element representation of 
the airfoil under consideration, the boundary conditions and data to specify 
additional loading conditions. A photograph of the airfoil looking down on 
the tip is shown in figure 1. The finite element representation of the air- 
foil is shown in figure 2 where both the nodal (grid) points and the elements 
are defined by numbers. The element numbers are in parenthesis to differen- 
tiate them from the numbers defining grid points. Note that the airfoil pro- 
jected plane has an aspect ratio (span/ chord) of approximately one. The 
finite element used is triangular, plane stress, anisotropic, and accomodates 
both bending and stretching. It is defined as element CTRIA2 in the NASTRAN 
L15.5 library of elements. The boundary conditions specified consist of 
fixing the three translations (displacements along x, y and z coordinate 
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directions) and the three rotations (slopes about the x, y and z coordinated 
directions) along the root (grid points 1 through 9, figure 2 ), The rotation 
about z is fixed at all grid points because of the relatively shallow camber 
of the airfoil, see figure 1. All other translations and rotations are free. 
Three additional inputs for loading conditions consist of: the airfoil rota- 
tional speed (250 revolutions per second, 2200 feet per second tip speed), 
the lamination temperature difference (-530^F) , and the combined temperature 
loading (operating temperature - lamination temperature). 

AIRFOIL INTERPOLATED DATA 

The input data described in the previous section Is operated on by the 
pre-processor (reference 2) to generate the following data: 

1. Coordinates of grid points 

2. Airfoil temperature at the grid points 

3. Airfoil pressure difference at the grid points 

4. Airfoil thickness at the grid points 

5* Ply number and stacking sequence at the grid points 

The data generated for items (1) to (5) above are in the form of internal 
tables (arrays) and are used as inputs by the pre -processor for subsequent 
operations. These data are also printed out in table form. Herein, the 
data generated for items (2) to (5) above are presented as contour plots in 
figures (3) to (6) in order to illustrate graphically their span-wise (root- 
tip) and chord-wise (leading edge-trailing edge) variations. Note that in 
figure 6 two numbers are given for each ply contour. The first number denotes 
the total number of plies while the second number in parentheses denotes the 
number of plies at 0®. Note, furthermore, the concentration of 0® plies in 
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the airfoil midchord region. Note also the similarity between ply contours 
(figure 6) and thickness contours (figure 5) , as is expected* 

AIRFOIL STIFFNESS AND THERMAL PROPERTIES 

The stiffness and thermal properties of the airfoil required for the 
finite element analysis constitute part of the airfoil structural response 
and are generated via the pre-processor (reference 2), The various airfoil 
stiffness and thermal properties generated along the x, y, z coordinates of 
the airfoil, figure 2, are the following: 

1. Moduli and Poisson’s ratio 

2. Thermal coefficients of expansion 

3* Membrane stiffnesses 

4* Bending stiffnesses 

5* Thermal forces 

6* Plane-stress stress-strain and temperature relationships 
compatible with the CTRIA2 NASTRAN finite element. 

The above properties are plotted in figures (7) to (14) to illustrate 
graphically their span-wise and chord-wise variations. The format in these 
figures is as follows. The properties are plotted along the chord at three 
span sections: root, midspan and tip. These are noted, respectively, by R, 
M, and T on the curves. The properties along the chord are plotted at the 
leading edge (noted L.E. on the curve), trailing edge (noted T.E. on the 
curve) and equidistant intermediate points (chord-fractions ) . 

Considerable insight on how the airfoil will respond when subjected to 
mechanical and/or thermal loads may be obtained by studying the properties 
plotted in figures (7) to (12). For example, in figure 7a, it is seen that 
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the modulus in the span direction is about twice as high at tnidchord than it 
is at the two edges. The converse is true for the in— plane shear modulus, 
figure 7c. Also the major Poisson ^s ratio is higher at the edges than it is 
at midchord and is large (greater than .60) when compared to homogeneous iso- 
tropic elastic materials (about .30). These large Poisson’s ratio values 
result from the dominance of the +45® plies at the edges. 

The thermal coefficient of expansion in the span direction, figure 8a, 
is negative at midchord and positive at the edges and is relatively small 
compared to structural metals. Because the coefficient is smaller in the 
midchord region than at the edges, increasing temperature will subject this 
airfoil to compression forces at the edges and tensile forces at midchord. 

The opposite effect will occur in the chord direction, figure Bb , since this 
coefficient of expansion is higher at midchord than it is at the edges. 

The membrane stiffness in the span direction > figure 9a, is 

greater at midchord than it is at the edges. This means that an applied load 
in the span direction, such as the centrifugal, will be carried mainly by the 
blade midchord region. This is so because there is a concentration of 0® 
plies in this region (figure 6). 

The bending stiffness in the span direction (D^^^) , figure 10a, shows 

similar variation as the membrane stiffness - ) . This will cause the 

'“11 

bending loads along the span to be carried primarily by the midchord portion 
of the airfoil. Note the variation of the bending stiffness in the chord 
direction (figure 10c) and the torsional stiffness (figure lOd) . Both of 
these are considerably lower near the edges than they are at the midchord 
region. Airfoils with stiffness distributions as shown in figures 10c and 
10 d are sxisceptlble to low frequency chordwise and torsional vibrations at 



- 8 - 


the edge regions and particularly near. the tip. 

Airfoils made from angleplied laminates will possess bending-twisting 
coupling. The distribution of stiffnesses which produce bending-twisting 
coupling are illustrated in figure 11. The magnitudes of the coupled stiff- 
nesses in figure 11 are considerably smaller than those in figure 10. Air- 
foils possessing coupling stiffnesses when subjected to bending loads will 
also twist. Airfoils with material coupling will vibrate in complex modes 
and not in the decoupled modes first bending, second bending, or first 
torsional as is generally considered to be the case with airfoils made from 
homogeneous isotropic materials. 

The distribution of the thermal forces is shown in figure 12. Note the 
sense of these forces is consistent with that mentioned in the discussion of 
the thermal coefficients of expansion, figure 8. The variation of the plane- 
stress stress-strain coefficients and the plane-stress thermal coefficients 
of expansion compatible with the CTRIA2 NASTRAN finite element are shown in 
figures 13 and 14 respectively. Note the variations of the thermal coeffi- 
cients in figures 14a and I4b are different from those in figure 8a and 8b 
because the former are evaluated at the element center while the latter at 
the node. Note also in figure 14c that the coefficient («^^) will produce 
shearing with changes in airfoil temperature. This shearing is caused by the 
presence of the in-plane stretching-shearing coupling coefficients (figures 
13e and 13f) . As a side note, the NASTRAN CTRIA2 plate element permits 
bending-twisting coupling via the moment-area integrals of the nonsyrametric 
terms and in the stress-strain relationships. 



AIRFOIL RESPONSE FROM VARIOUS LOADING CONDITIONS 


The airfoil under consideration was subjected to the various load 
conditions listed in Table 3* The response from these load conditions con- 
sists of: root reactions, displacements, composite surface stresses, through 
the thickness shear forces, maximum stressed ply stresses, interply delami- 
nation, and margins of safety* These responses were determined using the 
analysis capability described in reference 2 and are discussed in detail in 
the sections that follow* 

Root Reactions 

The force root reactions are presented graphically in figure 15. Note 
in this figure that the reaction forces due to centrifugal forces dominate 
those from the other loading conditions* The magnitudes of these reaction 
forces are used to design the airfoil root attachment. Note the high mag- 
nitude of F near the trailing edge, figure 15c. This indicates that forces 
z 

which tend to untwist the airfoil are set up* 

Corresponding plots for the root reaction moments are shown in figure 
16. The reaction moments due to centrifugal force dominate those from the 
other loading conditions. 

Airfoil Tip Displacements 

The airfoil tip displacements due to the various loading conditions 
examined are presented graphically in figure 17. Note the scale for the 
ordinates in this figure are irregular; however, the scale varies linearly 
between the minimum and maximum values. The displacement in the span 
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directlon indicates radial growth, in the chord direction, uncambering, and 
in the z-direction (thickness direction, approximately) out-of-plane bending 
and twisting* 

The results in figure 17 are useful in determining pre-twist, pre- 
camber and pre-tilt in order to keep the resulting displacements within 
design tolerance requirements. Note the tip displacements due to centrifugal 
load dominate those from the other loading conditions. Note also that the 
chord direction and ^bending” tip displacements due to pressure are in the 
opposite direction as compared with those due to centrifugal load. Further-- 
more, those due to thermal load are negligible compared to those from the 
other loading conditions. 

Composite Surface Stresses 

The composite surface stress in the span direction presented 

graphically in figures 18 (pressure surface) and 19 (suction surface). In 

these figures, the a variation due to various loading conditions is shown 
^ cxx 

along the chord at three sections along the span: root (R) , midspan (M) , and 
tip (T) . The following are observed from figures 18 and 19: 

1, The thermal stress is negligible compared to that from the 
centrifugal loading, 

2. The pressure stress is about 20 percent of the centrifugal stress. 

3- The Stress reaches its maximum magnitude (125 ksi) on the 

pressure surface at the root near the leading edge, figure 18d. 

The stress ^ minimum magnitude (60 ksi) at midspan in the 
midchord region (figure 19d) , 


4 , 
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The variation of the composite surface stress in the chord direction 
is presented graphically in figures 20 (pressure surface) and 21 
(suction surface)* The maximum magnitude of this stress (100 ksi) is on the 
pressure surface at the leading edge in the root> figure 20d. Everywhere 
else the magnitude is less than 20 ksi. The high magnitude of 100 ksi is 
also caused iri part by the specification of the boundary conditions. Allow- 
ing the airfoil to move along the root chord near the leading edge will 
decrease this magnitude substantially. 

The variation of the composite surface shear stress (a ) on the 

c xy 

pressure surface due to various loading conditions is presented graphically 
in figure 22. The corresponding variation on the suction surface is pre- 
sented in figure 23. The following are observed from the results in figures 
22 and 23: 


1. The shear stress due to centrifugal load is the dominant stress. 

2. The pressure surface shear stress (figure 22) has a maximum 
magnitude (60 ksi) in the leading edge near the root and a minimum 
magnitude (-40 ksi) in the leading edge at the midspan region. 

3. The suction surface shear stress (figure 23) has a maximum 
magnitude (30 ksi) in the leading edge at the root and a minimum 
magnitude (15 ksi) in the trailing edge at the midspan region. 

4. The shear stress has the same sign in the leading edge on both 
the pressure and suction surface. Similarly, the sign is the same 
in the trailing edge. This indicates that the airfoil is under- 
going untwisting. 
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Through-the-Thickness Shear Forces 

The variation of the shear force through the airfoil thickness on 

a section normal to the span due to various loading conditions is presented 

graphically in figure 24. As can be seen in this figure, the dominant shear 

force is that due to the centrifugal load. The combined shear force has a 

maximum magnitude of 2000 pounds/inch in the midchord region at the root, 

figure 24d. The maximum shear stress due to this force is given by the 

formula 1.5Q /t where t is the airfoil thickness at that section. Using 
xz 

2000 pounds /inch for and 0.5 inch for t yields a maximum shear stress of 
6000 psi. This value is about 40 percent of interlaminar shear strength as 
determined by the short beam shear test for the composite system used in the 
airfoil under consideration. 

The corresponding variation of the shear force through the thick- 

ness on a section normal to the chord is shown in figure 25. The maximum 
shear stress corresponding to the maxiTnum shear force, figure 25d, is 

7800 psi- Since both Q and Q cause stresses which are maximum at the 
^ xz yz 

same plane, they could cause delamination, similar to that observed in a 
short-beam-shear (interlaminar shear) test. 


MAXIMUM-STRESS ED-PLY STRESSES 

The most effective way to assess the adequacy of an airfoil laminate 
configuration with respect to stress is at the ply level. This is because 
the combination of +45® and 0® plies varies considerably throughout the 
airfoil and a large number of laminates need to be tested, at combined 
loads, if the assessment is to be made at the laminate level. 
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CoTnbined Loading at Operating Temperature 

For the airfoil under consideration the ply stresses of the maximuTn- 
s tressed— ply due to combined load (without residual stress) are sumniarized 
in Table 4. Data are con^jiled for several chord fractions near the root» 
midspan, and near the tip sections. Referlng to figure 2, these are 
respectively; nodes 10—18, 37—45, and 64—72. Table 4 contains the following 
information: 

1. Span section and chord-f raction first column 

2. Ntimber of plies at chord-fraction, second column 

3. Maximum stressed ply, third column 

4. Ply stress: along fiber * transverse to fiber (^^ 22 ^ 

intralaminar stress 

5. Ply combined-stress margin of safety (MDS) based on modified 
distortion energy principle refs. 3, 4, or 5, column 7 

6. The last three columns have information with respect to interply 
delaraination ; adjacent plies, angle change (A8) and margin of 
safety (MOS) (ref. 5), respectively. 

7. Unidirectional composite fracture stresses at the bottom of the 
table from which the ply MOS is calculated as discussed in 
refs- 3, 4, or 5. 

Note negative MOS is Tables 4 and 5 denote ply failure, and thus a critically 
stressed ply. The following are observed from the results in Table 4; 

1. The maximum stressed plies are within the shell; the shell 
includes the first 14 plies, and the last 14 plies at the 


maximum blade thickness. 
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2. Negative MOS*s are indicated in plies in the leading and trailing 
edges at the root and midspan and in the trailing edge at the tip. 
Comparing the individual ply stresses with the unidirectional 
composite fracture stresses it is seen that the negative MOS is 
caused primarily by either high intralaminar shear stress or high 
transverse stress in all cases except in the trailing edge at 
midspan. At this latter point the negative MOS is caused by high 
longitudinal compressive stress. 

3. The core (0*^) plies carry less load compared to shell (+45^) plies 
at the same chord fraction since the high stressed plies are in the 
shell. Though detail data are not presented here, the <^£11 
magnitudes in the core plies are less than 75 percent of those in 
the maximum-stressed shell plies and in many cases less than 50 
percent. Therefore the core plies are used inefficiently in this 
airfoil laminate configuratlon- 

4. Interply delamination or near delamination (MOS :tr0) is indicated 
in adjacent plies within the shell and in adjacent plies between 
shell and core. This means that the core-shell laminate configura- 
tion for airfoils in general may be susceptible to interply delami- 
nation . 

Lamination Residual Stress 

The ply stresses of the maximum-s tressed-ply due to lamination residual 
stress is induced by the temperature difference between cure and room 
temperatures (-530®F). The ply stresses from this source are summarized in 
Table 5. The following are observed from these results: 
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1. The maximum stressed plies are within the shell- 

2. Negative MOS's are indicated along the trailing edge- In all 
cases this is due primarily to a high transverse tensile stress. 

At the midspan and tip section, there is also a high longitudinal 
compressive stress. 

3. Interply delamination or near delamination (MOSCO) is indicated 
within the shell. This means that this laminate configuration is 
susceptible to transverse ply (transply) cracking and interply 
delamination at these locations after molding but without additional 
loading. 

Combined Loading with Residual Tliermal Stresses at Operational Conditions 

The temperature load condition for this case was obtained by using 
the operational temperature (figure 3) minus the lamination temperature 
(600®F). It is a combined case with the pressure and centrifugal forces 
also acting. This case, therefore, will include the ply residual stresses 
at the operational temperature and the stresses due to operating conditions, 
The ply stresses of the raaxiraum-s tressed-ply due to combined stresses are 
summarized in Table 6. The following are observed from these results: 

1, The maximum stressed plies are within the shell. 

2. Negative MOS’s are indicated along the leading edge at the root 
and midspan section and all along the trailing edge. The negative 
MOS is caused primarily by intralaminar shear or high transverse 
tensile stress except at the trailing edge of the midspan section. 

At this point, the negative MOS is caused by high longitudinal 


compressive stress. 
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Note that the difference between the data in Table 6 and that in Table 4 is 
the ply residual stresses at operating conditions. Up to this point the 
results presented and discussed are based on a particular airfoil laminate 
configuration. However the data generated can be used to modify the laminate 
configuration for more efficient cor^osite material utilization. 

AIRFOIL VIBRATION 

The vibrating of the airfoil under consideration was determined at 
zero-speed (no centrifugal force). The first 9 fundamental modes in cycles 
per second (CPS) are listed in ascending order in Table 7. As can be seen 
in this table the magnitude of the fundamental frequencies ranges from 234 
CPS to 2495 CPS. 

The reaction forces at the root due to first four fundamental frequen- 
cies are presented graphically in figure 26. As can be seen in figure 26, 
the reaction forces are of cyclic nature ranging from: -8000 to 7000 pounds 

for F , -3000 to 4500 pounds for F , and -1500 to 3000 pounds for F . These 

X y 2 

magnitudes are comparable to those of the root reaction forces due to centri- 
fugal load, figure 15. The vibration reaction forces are in the high cycle 
fatigue regime and could produce difficulties if the root attachment Is not 
designed to accomodate these cyclic loads. 

SUMMARY OF RESULTS AND CONCLUSIONS 

The major results and conclusions of an investigation to theoretically 
examine the structural response of composite airfoils for high-tip-speed 
compressor fan blades for aircraft turbine engines subjected to various 
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loading conditions (thermal, pressure, centrifugal, and combinations of 
these) are as follows: 

1. The root reactions, tip displacements and composite stresses due 
to thermal and pressure loading conditions are negligible compared 
to those produced by the centrifugal loading. 

2. The combined loading condition produced high stresses in the shell 

plies of the core-shell airfoil laminate configuration core, 

+45° shell) considered. These high stresses could produce ply 
failure in the shell plies in the leading and trailing edge regions 
at the root and inidspan sections, and in the trailing edge at the 
tip section. 

3. The core plies carry less load compared to shell plies resulting 
in inefficient use of the core plies in the core-shell type of 
configuration. 

4. The combined loading case could produce interply delaraination due 
to relative rotation in adjacent plies in the shell and between 
shell and core at the midchord and in the leading and trailing 
edge regions at the root and midspan sections. 

5. The combined loading case produced high through-t he- thickness 
shear stresses at the midchord in the root section. These stresses 
have sufficiently high magnitudes to produce interply delamination 
similar to that observed in short-beam-shear (interlaminar shear) 
tests . 

6. The root reactions due to various vibration modes are cyclic and 
are comparable to those of the combined load case. These reactions 
may subject the root attachment to high cycle fatigue. 
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7. The stiffness and therm^ properties of the airfoil section may be 
used to obtain an a priori assessment of the airfoil structural 
response. 

8. The core-shell type of laminate configuration for airfoils for 
high-tip-speed coiiqjressor fan blade applications appears to be 
inefficient with respect to high ply stresses near the edges and 
interply de lamination between shell and core. 

9. The lamination residual stresses in the plies at room temperature 
could produce ply failure (transply cracks) and interply delamina- 
tion in the shell plies along the trailing edge. These failures 
are caused primarily by high transverse and intralaminar shear 
stresses. 

10. The combined loading at operational temperature causes ply failure 
in the shell plies along both the leading and trailing edges. 

These failures are caused primarily by high transverse and intra- 
laminar shear stresses. 

11. The results presented herein may be used as a guide for high-tip- 
speed airfoil designs using composites with respect to: stiffness, 
thermal loads, root reactions, composite and ply stresses, tip 
displacements and fundamental frequencies. 

12. Preliminary designs for high-tip-speed composite airfoils for suf- 
ficient strength and minimum tip displacements may be based only on 
the centrifugal load since the structural response from this load 
dominates all the other loading conditions. 
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NOMENCLATURE 


A axial stiffness array -position denoted by numerical subscription 

c 

C,^ subscript - composite property 

bending stiffness array - position denoted by numerical subscripts 

F reaction force - direction denoted by subscript; terminal point 

of ply-type 

FVR fiber volume ratio 

G array of anisotropic stress-strain relationships for the NASTRAN 

plate elements - position denoted by numerical subscripts 

I starting point of ply-type 

L.E. leading edge of airfoil 

subscript - ply property 

M midspan section; moment -direction denoted by numerical subscripts 

MOS margin of safety 

P pressure 

Q shear force through airfoil thickness -surface and direction denoted 

by subscripts 

R root section 

T temperature 

T.E. trailing edge of airfoil 

t airfoil thickness 

VVR void volume ratio 

x,y,z airfoil structural axes, x-along span; y-along chord; z-along thick 
ness, subscripts — corresponding coordinate directions 

1,2,3 subscripts - array position or ply principal direction with 
along fiber 

ZL Z - coordinate to pressure surface 

ZU Z - coordinate to suction surface 

a thermal coefficient of expans ion- direct ion denoted by numerical 

subscript 
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AG interply relative rotation 

0 stress-type and direction denoted by subscripts 



OJLOUJUJUJUOLiJUiOJiOKJrOhONJNJtO rotOMM 


TABLE 1 INPUT DATA 


BLADE GEOMETRY, TEMPERATURE AND PRESSURE 


SECTION 

1 

1 

1 

1 

1 

1 


X 

Y 

ZU 

ZL 

T 

DELTA-P 

RADIUS 

CHORD 

SUCTION SURFACE 

PRESSURE SURFACE 

TEMPERATURE 

PRESSURE DIFFERENCE 

9.38 

.000 

• 060 

.000 

181. 

10.5 

9.38 

1.00 

.500 

.200 

175. 

10.0 

9.38 

2.00 

.800 

.400 

190. 

9.00 

9.38 

3.00 

.955 

.440 

205. 

7.00 

9.38 

4.00 

.910 

.404 

212. 

5.50 

9.38 

5.00 

.770 

.310 

220. 

4.50 

9.38 

6.00 

.350 

.121 

227. 

4.00 

9.38 

6.60 

.062 

.000 

235, 

3-50 

13.0 

-.150 

-.900 

-.950 

285. 

14.5 

13.0 

.000 

-.850 

-.900 

275. 

14.5 

13.0 

1.00 

-.400 

-.550 

265. 

12.5 

13.0 

2.00 

• 050 

-.200 

275. 

9.50 

13.0 

3.00 

.650 

.200 

285. 

7.50 

13.0 

4.00 

1.05 

.590 

287. 

6.00 

13.0 

5.00 

1.35 

.980 

289. 

6.00 

13.0 

6.00 

1.55 

1.40 

292. 

7.00 

13.0 

6.65 

1.75 

1.70 

296. 

7.20 

16.4 

-.300 

-1.31 

-1.35 

443. 

15.0 

16.4 

.000 

-1.20 

-1.25 

415. 

14.0 

16.4 

1.00 

-.700 

-.773 

386. 

13.0 

16.4 

2.00 

-.150 

-.350 

386. 

10.0 

16.4 

3.00 

.450 

• 050 

386. 

7.00 

16.4 

4.00 

1.00 

.650 

386. 

8.00 

16.4 

5.00 

1.50 

1.30 

386. 

12.0 

16.4 

6.00 

2.17 

2.07 

386. 

14.0 

16.4 

6.50 

2.50 

2.45 

388. 

12.0 


NOTE ; Dimensions are as follows : 
Distance Inches 

Temperature ‘^F 
Pressure PSI 





TABLE 2 INPUT DATA OF COMPOSITE TYPE AND PLY GEOMETRY 


THICKNESS 

PLY 

DESIG- 

X-COORDINATE 

Y-COORDINATE 

TYPE OF COMP. 

PLY 

PLY 

PLY 

%F 

ORIENT 

NATION 

%i 

%F 

%i 

%F 

FIBER/MATRIX 

THICKNESS 

VVR 

FVR 

0*0 20.0 

45.0 

1.0 

.0 

100.1 

.0 

100.1 

MOD2 

POLY 

.007 

.000 

.570 

20.0 50.1 

.0 

2.0 

.0 

100.1 

.0 

100.1 

M0D2 

POLY 

.007 

.000 

.570 


%I 

Starting point 

of ply-type 

%F 

Terminal point 

of ply-t3^e 

WR 

Void Volume Ratio 

FVR 

Fiber Volume 

Ratio 


TABLE 3 SUMMARY OF CASES EXAMINED 


CASE 

1 

2 

3 

4 

5 

6 


DESCRIPTION 
THERMAL LOAD ONLY 
PRESSURE LOAD ONLY 
CENTRIFUGAL FORCE ONLY 
COMBINED . CASES 1-3 
RESIDUAL STRESS 
OPERATING THERMAL STRESS 


7 


VIBRATION MODES 



TABLE 4 MAXIMUM-STRESSED- PLY STRESSES AT SELECTED CROSS SECTIONS DUE TO COMBINED LOADING 




PLIES 

MAX. 


PLY 

STRESSES (KSI) 


PLY 

RELATIVE 

ROTATION 

CHORD 


PER 

STRESSED 





COMB, STRESS 

ADJ. 

AO 

COMB.STRES 

FRACTION 


NODE 

PLY 

0)111 

(jZ22 

ail2 

MOS 

PLIES 

DEG 

MOS 

ROOT 

L.E. 

7 

7 

- 96.4 


30.1 

8.2 

- 23. 

6/7 

-.025 

-2.93 

SECTION: 

1/8 

30 

30 

15.8 

- 

3.1 

16.2 

.081 

16/17 

,n052 

,190 


1/4 

45 

45 

65.3 

- 

1.4 

12.4 

.299 

31/32 

.0051 

,206 


3/8 

58 

58 

114. 

- 

.2 

9.4 

.246 

51/52 

-.0064 

.012 


1/2 

71 

71 

103. 

- 

0 

3.4 

,592 

70/71 

-.0059 

.092 


5/8 

70 

2 

56.4 

- 

1.3 

- 1.6 

.845 

1/2 

.0046 

.295 


3/4 

64 

63 

31.4 

- 

2.1 

.9 

.924 

62/63 

.0039 

.395 


7/8 

37 

15 

21.2 

- 

4.0 

.1 

.904 

35/36 

.0034 

.467 


T.E. 

8 

1 

31.4 

- 

8.2 

- 36.3 

- 3.66 

2/3 

-.0067 

- .039 

MID - 

L.E. 

6 

1 

-159.3 


5.3 

- 22.6 

- 1.52 

1/2 

.0086 

- .331 

SPAN : 

SECTION 

1/8 

16 

1 

- 51.6 

- 

3.3 

- 8.9 

.667 

1/2 

.0052 

.204 

1/4 

27 j 

27 1 

78.1 

- 

1.0 

3.1 

.723 

26/27 

-.0054 

.163 


3/8 

42 1 

42 

108.2 

- 

.9 

1.6 

.545 

35/36 

-.0063 

.024 


1/2 

60 

60 i 

87.1 

- 

.7 

0 

.711 

59/60 

-.0056 

.131 


5/8 

61 

47 1 

67.1 

- 

2.7 

.5 

.751 

59/60 

.0042 

.356 


3/4 

51 

1 

56.5 

- 

4.1 

.2 

.749 

2/3 

-.0042 

.343 


7/8 

26 

25 1 

34.0 

- 

2.0 

8.5 

.680 

2/3 

-.0050 

.235 


T.E. 

6 

1 

■ - - ‘ 

2 ! 

-241.8 

- 

6.5 

7.1 

- 1.21 

2/3 

-.0095 

— .465 





TABLE 4 CONTINUED 


CHORD 

FRACTION 

PLIES 

PER 

NODE 

MAX. 

STRESSED 

PLY 

PLY STRESSES (KSI) 

PLY RELATIVE 

ROTATION 

aall 

aZ22 

0^12 

CO MB. STRESS 
MOS 

ALJ. 

PLIES 

AO 

DEG 

COMB.STRES 

MOS 

TIP ^ L,E, 

5 

2 

24.9 

- 4.6 

3.0 

.842 

1/2 

.0062 

.049 

SECTION' 

8 

8 

- 2,3 

- 4.8 

4.6 

.871 

1/2 

.0049 

.250 

1/4 

17 

17 

2.7 

- 4.3 

- 5.5 

.845 

16/17 

-.0046 

.296 

3/8 

33 

19 

32.8 

- 5.2 

“ .4 

.819 

32/33 

-.0049 

.247 

1/2 

51 

37 

16.4 

- 4.8 

• ^ ! 

.894 

50/51 

-.0051 

.215 

5/8 

48 

34 

25.8 

- 4.8 

0 

.863 

47/48 

-.0040 

.383 

3/4 

33 

19 

36.0 

- 5.2 

.3 

.802 

32/33 

-.0049 

.245 

7/8 

17 

17 

32.5 

- 3.8 

- 8.8 

.614 

16/17 

-.0058 

.110 

T.E. 

5 

5 

85.0 

.4 

- 34.4 

- 3.181 

4/5 

-.0053 

.188 


NOTES: Unidirectional Composite 

Material Fracture Stresses 


Longitudinal Tensile 

170 

KSI 

Longitudinal Compression 

150 

KSI 

Transverse Tensile 

9 

KSI 

Transverse Compression 

20 

KSI 

Intralaminar Shear 

17 

KSI 

Interply Delamination 

Critical Angle A 8 ^ 

.0065 DEG 
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TABLE 5 MAXIMUM -STRESS ED- PLY STRESSES AT SELECTED CROSS SECTIONS DUE TO LAMINATION RESIDUAL STRESS 


CHORD 

FRACTION 

PLIES 

PER 

NODE 

MAX. 

STRESSED 

PLY 

PLY STRESSES OKSn 

PLY 

RELATIVE ROTATION 

Oj^ll 


aU2 

COMB.STRES 

MOS 

3 ADJ. 
PLIES 

AO 

DEG 

COMB. S TREE 
MOS 

ROOT 

L.E, 

7 

6 


1.37 

6.06 

- 2.61 

.578 

1/2 

-.0062 

.034 

SECTION: 

1/8 

30 

16 

- 

8.61 

5.98 

.02 

.552 

1/2 

-.0059 

.082 


1/4 

45 

31 


4.95 

6.04 

.00 

,568 

43/44 

".0059 

.083 


3/8 

58 

44 

- 

3.81 

6.08 

.01 

.568 

2/3 

.0059 

.073 


1/2 

71 

15 

- 

4.15 

5.94 

.02 

.586 

69/70 

-.0059 

.079 


5/8 

70 

15 

- 

5.52 

6.11 

.11 

.555 

68/69 

-.0061 

.061 


3/4 

64 

15 


5.75 

6.08 

.04 

.558 

62/63 

-.0059 

.077 


7/8 

37 

15 

- 

7.44 

5.99 

.00 

.559 

2/3 

.0059 

.084 


T.E. 

8 

8 

1 

- 

36.0 

227.8 

.98 

-580.0 

7/8 

.2357 

-35.4 

MID - 

L.E. 

6 

1 

_ 

4.14 

5.66 

.00 

.623 

2/3 

.0060 

.067 

SPAN 

SECTION: 

1/8 

16 

1 

- 

3.06 

5.63 

i 

.04 

.633 

2/3 

.0059 

.082 


1/4 

27 

27 

- 

2.65 

5.64 

.08 

.634 

26/27 

.0059 

.084 


3/8 

42 

15 

- 

3.18 

5.66 

.00 

.628 

41/42 

.0059 

.082 


1/2 

60 

15 

- 

3.04 

5.67 

.00 

.628 

1/2 

-.0059 

,084 


5/8 

61 

47 

- 

3.29 

5,66 

.00 

.627 

1/2 

-.0059 

.082 


3/4 

51 

1 2 

1 

- 

3.00 

5.62 

.07 

.634 

1/2 

-.0059 

.084 


7/8 

26 

1 

! 2 

- 

3.02 

5.63 

.25 

.633 

1/2 

-.0059 

.082 

1 

T.E, 

6 

6 

"259.7 

504.8 

1.88 

-2900 

5/6 

.5325 

-81.21 
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TABLE 5 CONTINUED 


CHORD 

FRACTION 

PLIES 

PER 

NODE 

MAX. 

STRESSED 

PLY 

PLY STRESSES (KSI) 

PLY 

RELATIVE 

ROTATION 

oJ.ll 

a£22 

all! 

COMB. STRESS 
MOS 

ADJ. 

PLIES 

A0 

DEG 

COMB.STRE 

MOS 

TIP L.E. 

5 

1 

- 2.60 

5.63 

- .09 

.636 

2/3 

.0059 

.086 

SECTION: 

8 

2 

- 2.68 

5.63 

.02 

.635 

1/2 

-.0059 

.084 

1/4 

17 

1 

- 2.75 

5.62 

.03 

.636 

2/3 

.0059 

.086 

3/8 

33 

15 

- 3.41 

5.63 

.01 

.630 

2/3 

.0059 

.084 

1/2 

51 

15 

- 3.05 

5.63 

.00 

.632 

2/3 

.0059 

.086 

5/8 

48 

15 

- 2.91 

5,61 

.00 

.635 

1/2 

-.0059 

.086 

3/4 

33 

1 

- 3.07 

5.61 

- .03 

.635 

2/3 

.0059 

.084 

7/8 

17 

1 

- 2.99 

5.63 

- .08 

.634 

2/3 

.0059 

,083 

T.E. 

5 

5 

-101.1 

218.5 

- ,13 

-550. 

4/5 

.2299 

-34.49 


NOTES: Unidirectional Composite 

Material Fracture Stresses 


Longitudinal Tensile 170 KSI 
Longitudinal Compression 150 KSI 
Transverse Tensile 9 KSI 
Transverse Compression 20 KSI 
Intralaminar Shear 17 KSI 
Interply De lamination 


Critical Angle A 0 .0065 DEG 
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TABLE 6 MAXIMUM- STRESS ED' PLY STRESSES AT SELECTED CROSS-SECTIONS DUE TO COMBINED THERMAL LOADING 
(THERMAL LOADING = OPERATING TEMP. - LAMINATION TEMP.) 


CHORD 

FRACTION 

PLIES 

PER 

NODE 

MAX. 

PLY STRESSES (KSI) i 

PLY RELATIVE ROTATION 

STRESSED 

PLY 

aJlll 

a%22 

oi!,12 

COMB . STRES 
MOS 

S ADJ. 

PLIES 

A0 

DEG 

COMB.STR2 

MOS 

ROOT 

7 

7 

- 89.6 

- 24.1 

85.4 

-24.0 

2/3 

.0013 

.787 

SECTION: 

30 

16 

- 67.0 

9.12 

.91 

- ,702 

28/29 

-.0057 

.120 

1/4 

45 

31 

- 61.9 

10.2 

1.61 

- .929 

41/42 

-.0064 

.009 

3/8 

58 

44 

- 32.6 

10.9 

2.56 

- .721 

•13/14 

-.0053 

.177 

1/2 

71 

57 

17.2 

8.57 

1.97 

.285 

61/62 

-.0064 

.015 

5/8 

70 

15 

13.2 

6.11 

- 1.23 

.642 

2/3 

.0052 

.195 

3/4 

64 

64 

6.77 

4.83 

- .60 

.766 

63/64 

.0046 

.294 

7/8 

37 

37 

- 8.69 

3.49 

- .62 

.832 

36/37 

.0041 

.372 

T.E. 

8 

1 

17.5 

- 3.25 

-31.2 

2.27 

3/4 

-.0053 

.183 

MID • L.E. 

6 

1 

-163.4 

.31 

-22.6 

- 1.9-7 

5/6 

.0037 

.429 

SPAN 

SECTION: ' 

16 

1 

- 54.6 

2.36 

- 8.93 

.417 

2/3 

,0050 

.225 

1/4 

27 

1 

- 40.1 

2.64 

- 7.75 

.551 

25/26 

-.0052 

.204 

3/8 

42 

41 

20.6 

7.57 

- 1.45 

.473 

38/39 

-.0064 

.017 

1/2 

60 

59 

84.3 

4.94 

.00 

.636 

58/59 

-.0052 

.202 

5/8 

61 

61 

28.3 

5.22 

- 3.60 

.748 

60/61 

.0038 

.420 

3/4 

51 

51 

13.0 

4.12 

- 4.25 

.791 

1/2 

-.0038 

.408 

7/8 

26 

26 

15.5 

4.35 

- 8.75 

.585 

1/2 

-.0045 

.297 

T.E. 

6 

2 

-250 . 

.88 

1 9.06 

J - 

- 2,06 

2/3 

-.0036 

.435 
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TABLE 6 CONTINUED 


CHORD 

FRACTION 

PLIES 

PER 

NODE 

MAX. 

STRESSED 

PLY 

PL^ 

1 STRESSES 

(KSI) 


PLY 

RELATIVE ROTATION 

ofcll 

a£22 

a£12 

COMB,STRES 

MOS 

s" ADJ. 
PLIES 

AO 

DEG 

COMB.STRE5 

MOS 

TIP ^ L.E, 

5 

1 

- 20,2 

2.97 

- 3.08 

.796 

2/3 

-.0022 

.654 

SECTION* 

8 

1 

- 2.53 

2.90 

- 3.29 

.864 

2/3 

.0026 

.604 

1/4 

17 

16 

- 12.2 

1.79 

5.25 

.846 

15/16 

-.0025 

.609 

3/8 

33 

32 

8.80 

1.52 

- 3.72 

.879 

19/20 

-.0028 

.560 

1/2 

51 

50 

- 16.0 

2.04 

1.64 

.903 

49/50 

-.0030 

.537 

5/8 

48 

48 

1.01 

1.92 

- 2.81 

.935 

34/35 

-.0022 

.660 

3/4 

33 

32 

- 5.39 

2.40 

5.68 

.816 

19/20 

-.0030 

.528 

7/8 

17 

16 

- 7.46 

2.98 

8.40 

.646 

15/16 

-.0035 

.463 

T.E. 

5 

5 

82.0 

5.24 

-34.6 

- 3.13 

3/4 

-.0028 

.570 


NOTES : 


Unidirectional Composite 



Material Practure Stresses 



Longitudinal Tensile 

170 

KSI 

Longitudinal Compression 

150 

KSI 

Transverse Tensile 

9 

KSI 

Transverse Compression 

20 

KSI 

Intralaminar Shear 

17 

KSI 

Interply Delamination 




Critical Angle A 0 ,0065 DEG 
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TABLE 7 FUNDAMENTAL FREQUENCIES OF AIRFOIL AT ZERO SPEED 


MODE 

1 

2 

3 

4 

5 

6 

7 

8 


FREQUENCY CYCLES/SECOND 
234 
736 
903 
1326 
1526 
1827 
1998 
2350 


9 


2495 



-8i47 


pa 


7 . 


5 . 


1 . 


0 

2-COORDll 
1+ OUTWAF 


Y-COORDINATE (IN. », CHORD DIRECTION 
















MODULUS (PSD 


L E. = LEADING EDGE 
T.E. - TRAILING EDGE 
R = ROOT SECTION 
M - MIDSPAN SECTION 
T-TIP SECTION 


15x10® 



DIRECTION. 


4x10® 



DIRECTION. 




CHORD FRACTION 


fc) SHEAR MODULUS (d) POISSON'S RATIO 

(x-y, FIG, 2). x-y, FIG, 21, 


Figure 7. - Airfoil moduli and Poisson's ratio as deter- 
mined by the pre-processor. 


LL = LEADING EDGE 
T. E. = TRAILING EDGE 
R - ROOT SECTION 
M = MIDSPAN SECTION 
T-TIP SECTION 


.50x10"® 




L E. - LEADING EDGE 
T. L - TRAILING EDGE 
R - ROOT SECTION 
M - MIDSPAN SECTION 
T • TIP SECTION 





CHORD FRACTION 


(bl THERMAL COEFFICIENT - CHORD DIRECTION. 

Figure 8. - Airfoil thermal coefficients of 
expansion as determined by the pre- 
processor. 


(cl CHORD DIRECTTON Id) IN-PLANE SHEAR 

Wc22*- (A(-33). 

Figure 9, - Airfoil membrane stiffness as determined by 
the pre- processor. 


BENDING STIFFNESS IN 10^ LB-IN. 


LE. - LEADING EDGE 
T.L - TRAILING EDGE 
R = ROOT SECTION 
M = MIDSPAN SECTION 
T - TIP SECTION 




IDcii). EFFECT 




CHORD FRACTION 

{cl CHORD DIRECTION (d) TORSIONAL 

(Dc22>. 


LE. - LEADING EDGE 
L E. = TRAILING EDGE 
R = ROOT SECTION 
M = MIDSPAN SECTION 
T - TIP SECTION 



CHORD FRACTION 

la) SPAN DIRECTION- (b) CHORD DIRECTION- 
BENDING-TWISTING BENDING-TWISTING 

COUPLING (Dqj 3 ). COUPLING ( 0 ^ 23 !* 

Figure 11. - Airfoil coupled bending -twisting stiffness 
as determined by the pre-prxessor. 


Figure 10. - Airfoil bending stiffnesses determined by 
the pre- processor. 


THERMAL FORCE LB/iN. 


LE. = LEADING EDGE 
1 . 1 = TRAILING EDGE 
R = ROOT SECTION 
M = MIDSPAN SECTION 
T = TIP SECTION 




LE. = LEADING EDGE 
T.E. - TRAILING EDGE 
R = ROOT SECTION 
M = MIDSPAN SECTION 
T-TIP SECTION 





CHORD FRACTION 

(c) CHORD DIRECTION (d) IN-PLANE SHEAR (G«). 
IG2?). 

Figure 13. - Plane-stress^ stress-strain relationships 
compatible with CTRIA2NASTRAN finite element 


Figure 12. - Airfoii thermai forces as deter- 
mined by pre-processor due to airfoii tem- 
perature shown in figure 2. 


STIFFNESS COEFFICIENT, # PSI 


5 000 1 — 


2000 ^ 



(e) COUPLES SPAN 
RECTfON STRETCHING 
WITH IN-PLANE 
SHEAR (Gi 3 ). 

1 . 25 1 — 



L.E. 1/4 1/2 3/4 T.E. 

CHORD FRACTION 


(f) COUPLES CHORD DI- 
RECTION STRETCHING 
WITH IN- PLANE 
SHEAR IG 23 L 

Figure 13, - Continued, 


L. E. = LEADING EDGE 
T.E = TRAILING EDGE 
R = ROOT SECTION 
M = Ml OS PAN SECTION 
T - TIP SECTION 



CHORD FRACTION 
ic) IN'PLANE SHEARING 
(012). 


Figure 14. - Plane stress thermal coefficients of 
expansion compatible with the CTRIA2 NASTRAN 
finite element. 




o 

o 


<c 

UJ 

a: 

8 


-20000^ 


J ^ ^ I 

LE 1/4 1/2 3/4 IE. 


-8000 I 


L 


L E. 1/4 1/2 3/4 


CHORD FRACTION 


(a) F^ 


J 

T. E. 


(b) Fy. 
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Figure 15. - Root reaction forces vs. blade chord for vari- 
ous load conditions (x,y,z airfoil coordinates, fig. 1 ). 
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Figure 17. - Airfoil tip displacements due to various loading conditions. 


Figure 16, - Root reaction moments vs, blade chord 
for various load conditions (x,y,z airfoil coordi- 
nates, fig. 1). 
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Figure 18. - Span direction stress various 

loading conditions. Pressure surface. 
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Figure 19. - Span direction stress locxx^ vari- 
ous loading conditions. Suction surface. 
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Figure 20. - Chord direction stress (oQyy) due to 
various loading conditions. Pressure surface. 
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Figure 2h - Chord direction stress (o^yy) due to 
various loading conditions. Suction surface. 
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Figure 2t - Shear stress to™) due to various loading 
cond i t j ons . P ress u re s u rrace. 
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Figure 23. ^ Shear stress (o™) due to various loading 
conditions. Suction surface^ 
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Figure 24. - Through -the-thickness shear force due 
to various loading conditions. 
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Figure 25, " Through -the-thickness shear force Oy^ 
due to various loading conditions. 
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Figure 2Cl - Root reaction force vs. blade chord for various 
vibration modes 1. 2. 3, 4 (x. y. z. airfoil coordinates 
fig. 1). 
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